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ABSTRACT 

The rapid time variabihty of gamma-ray bursts imphes the sources are very compact, and the 
peak luminosities are so high that some matter must be ejected at ultra-relativistic speeds. The 
very large Lorentz factors of the bulk flow are also indicated by the very broad and hard spectra. 
It is natural to expect that when the relativistic ejecta interact with the interstellar matter a 
strong synchrotron radio emission is generated, as is the case with supernova remnants and radio 
galaxies. We estimate that the strongest gamma-ray bursts may be followed by radio transients 
with peak fluxes as high as 0.1 Jy. The time of peak radio emission depends on the distance scale; 
it is less than a minute if the bursts are in the galactic halo, and about a week if the bursts are at 
cosmological distances. 

Subject headings: cosmology: theory — gamma rays: bursts — radio continuum: general — relativity 
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1. INTRODUCTION 

The isotropic distribution of gamma-ray bursts over the sky and non-uniform distribution in 
distance (Meegan et al. 1992) indicates that we are located at or near the center of a spherical 
and bound distribution of sources. This is easiest to understand if the sources are at cosmological 
distances (Paczyhski 1991, Dermer 1992, Mao & Paczynski 1992, Piran 1992). The inferred peak 
luminosity is ~ 10^^ erg s"^ > 10^^ Lq, which is in excess of Eddington luminosity for any 
object with M < 10^^ Mq, and is highly super-Eddington for any object that might generate a 
variability on a sub-millisecond time scale. This implies that the radiation pressure must drive 
an ultra-relativistic wind (Paczyhski 1990, and references therein). This conclusion remains valid 
even if the radiation is beamed into a small solid angle. The observed spectra are very broad and 
non-thermal, extending to at least 10 MeV and in some cases beyond 100 MeV and showing no 
evidence for a pair creation cut-oflF at ~ 511 keV (Schaefer et al. 1992, and references therein). This 
is independent evidence for ultra-relativistic outflow (Goodman 1986, Paczyhski 1986, Fenimore 
et al. 1992, 1993). 

There are two other types of gamma-ray sources with spectra as broad as burst spectra: 
some radio pulsars (Ruderman &; Cheng 1988) and some AGNs (Hartman et al 1992, Dermer & 
Schlickeiser 1992.) Both classes of object are also strong radio-emitters. There may be a very 
general reason for this, as both radio and gamma-ray emission require the presence of relativistic 
particles. Therefore, it is natural to expect that gamma-ray bursts may give rise to radio emission 
as well (Paczyhski 1992). 

There is plenty of evidence that all explosive events in astrophysics give rise to synchrotron 
radio-emission when the ejecta interact with the ambient medium. Radio jets generated by active 
galactic nuclei create 'hot spots' when they are stopped by extragalactic matter, and the products of 
the collision create giant radio-lobes, which outshine the core sources at low frequences. Supernova 
remnants are also powerful radio-sources. 

The most spectacular radio supernova known to date is SN 1986 J in NGC 891, at a distance 
of ~ 12 Mpc. The peak intensity was ~ 100 mJy over the frequency range 1-10 GHz (Rupen 
et al 1987) and the emission lasted for a few years. The total radio energy generated by the 
supernova was Eradio ~ 10"^^' erg, and a surface brightness of Tb ~ 5 x 10^^ K was reached, close to 
the Compton limit of 10^^ K (Kellerman & Pauliny-Toth 1969). The kinetic energy of supernova 
ejecta is typically Ekin ~ 10^^ erg, hence the radio efficiency of SN 1986J was ~ 10~^. We note 
that a ten minute observation with the Very Large Array (VLA) reaches sensitivity limits ~ (0.08, 
0.1, 0.06, 0.2, 0.3) mJy at frequencies (1.4, 5.0, 8.4, 15, 23) GHz, respectively (Very Large Array 
Observational Status Summary, 26 May 1993, page 9, NRAO Publication). Therefore, a supernova 
like SN 1986J is detectable out to a distance d « 0.4 Gpc. 

A 'typical' gamma-ray burst radiates ~ 10^^ ergs in the 20 keV - 2,000 keV range, assuming 
spherical symmetry. Adopting ~ 10% gamma ray efficiency, the total energy in the burst may be 

as large as Eq ^ lO^"' ergs. If the ejecta of some bursts are as efficient radio emitters as SN 1986J 
then they may be detectable with the VLA out to 1 Gpc and beyond. 

In the following section we make the case for radio emission from gamma-ray bursters some- 
what more quantitative, and in the last section we discuss the observable consequences of our 
scenario. 



2. RADIO FIREBALL 
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There are ~ 40 gamma ray bursts per year per 4tt steradians out to a redshift z ~ 0.2 (Mao 
& Paczynski 1992). Every year the positions for a handful of the strongest bursts are determined 
with ~ 1 arc minute accuracy using the interplanetary network (IPN) (Cline et al. 1992, and 
references therein). A typical distance to these bursts is hkely to be 0.5 Gpc or so. 

The following sequence of physical conditions seems to be fairly general, and does not depend 

on the mechanism responsible for the primary energy source or on the nature of the bursting object. 
The energy density at the source is so enormous that the optical depth to all elementary processes is 
very high and the conditions are close to local thermodynamic equilibrium, which implies that the 
gamma-ray spectrum in the region is close to blackbody form (Goodman 1986, Paczynski 1986). 
However, the observed spectra are very broad and very different from any blackbody (Schaefer et 
al. 1992), and therefore the observed emission cannot be due to the original fireball. Fortunately, 
even a small 'baryon loading' is sufficient to convert almost all the original energy into kinetic 
energy of the relativistic wind (Paczynski 1990, Shemi & Piran 1990). In the observer's frame the 
fireball looks like a thin spherical shell expanding ultra-relativistically (Blandford &: McKee 1976, 
Shemi &: Piran 1990). The shell cools rapidly because of adiabatic expansion. 

The absence of the pair creation cut-off in the observed spectra implies that the bulk Lorentz 
factor at the time of gamma-ray emission is very large, Fj > 10^ (Fcnimorc ct al 1992, 1993). The 
initial bulk Lorentz factor of the fireball Fq can only by larger. The collision between the ejecta 
and circum-source matter, and later interstellar and/or intergalactic matter, will gradually lower 
the bulk Lorentz factor, while converting the bulk kinetic energy into random energy of relativistic 
particles as described by Rees and Meszaros (1992) and Meszaros and Rees (1993) in their model 
of gamma-ray bursts. This 'randomized' energy is partly radiated away and partly used up for 
further expansion, i.e. it is re-converted into kinetic energy of the bulk fiow. Ultimately, after a 
long enough time interval, all the initial energy is either radiated away or 'wasted' in the Hubble 
expansion of the universe. 

Assuming that the emission is incoherent, one can estimate the synchrotron and inverse Comp- 
ton radiation from the fireball (van der Laan 1966, Rees 1967, Gould 1979). There is an upper 
limit Tb ^ 10^^ K to the brightness temperature of an incoherent synchrotron source (Kellerman 
& Pauliny-Toth 1969, Kellerman 1974). In the observer's frame this limit is increased by the 
Lorentz factor of the bulk fiow. Beyond this limit the energy density in radiation {urad) exceeds 
that in magnetic fields (Ug ), resulting in the so-called "Compton catastrophe," where the emitting 
electrons lose energy very rapidly through inverse Compton scattering. Such conditions are likely 
to occur early in the evolution of the blast wave, while all energy densities are still very high. 
If they do occur, the resulting inverse Compton emission may produce the observed gamma-ray 
bursts - a scenario similar to the model of Rees and Meszaros. When the expansion continues 
and the energy density falls below some critical level, synchrotron emission dominates over inverse 
Compton losses. This is the regime we are interested in. 

The following analysis will be done assuming spherical symmetry of the fireball. We begin with 
the injection of energy Eq into a small volume (rg <^ 300km) with the initial rest mass Mq. The 
initial fireball quickly evolves into a thin shell expanding with bulk Lorentz factor Fq = Eq/Mqc^ . 
By the time it has expanded to ^ Fgro, it becomes cold in the comoving frame, i.e. the internal 
Lorentz factor 7^ ~ 1. We assume that the ambient medium is at rest and has a uniform mass 
density p. As the blast wave expands, this material is swept up and mixes with the shell, which is 
reheated and decelerated in the process. Parametrizing the swept up rest mass by M = Mo(l -I-/), 
we can write energy and momentum conservation as 



ro + / = (1 + fhiV, 
ro/3o = (1 + /)7ir/3, 



(1) 

(2) 
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where F = (1 — is the bulk Lorentz factor, (3 = v/c as usual, and 7^ is the internal Lorentz 

factor of the particles moving randomly within the expanding shell. Solving these gives F and 7j 
as functions of /: 

To + f ^Vl + 2rof + f^ /?. = (l-7.-Y/^ (4) 

(i + f)r 1 + f ' A u 7. j • w 

There are three diflFerent measures of time which are relevant. The first, t, is measured in the 
rest frame of the burster, so that 

r^ct, M = Mo(l + /) = Mo + Y(ct)3p, (5) 

where r is the radius of the shell. The second, tco, is measured in the frame comoving with the 
shell, so that 

(6) 

The third, t^, is measured in the terrestrial observer's frame and must account for apparent 
superluminal motion, so that 

tQ)=t -^cosO. (7) 

We may write these in terms of / as 

too ~ (1 + o.32ro/)'/' , o</^ro (9) 

te^4^a + 0-5^of) . 0</^ro (10) 

Here the approximate form for tco interpolates between results accurate for f <^ T^^ and for 
Tq^ <C / <C Fq, while tq) is calculated for 6 = 0, i.e. along the line of sight. It is convenient to 
develop our model first in terms of / and later transform to . 

The expanding shell as seen in its comoving frame will expand at a rate comparable to the 
sound speed c/\/3, so that the thickness is Ar « ctco- The shell's surface is perpendicular to the 
direction of expansion, so it is the same in the comoving and stationary frames, A = 47r(ct)^. Thus, 
in the frame co-moving with the shell its volume is given as 

V^47rict)\ctco)^^f{l + 0.32roff' , (11) 

pro 

and its total energy as 

Ei = UiV = Mo(l + /)7zc' = Moc^ (l + 2ro/ + ff^ , (12) 
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where Ui is the total energy density. 

The expanding sheh is reheated to relativistic energies and the internal Lorentz factor may 
be as large as ji,max ~ (ro/2)^^^ (cf. eq. 4). These are ideal conditions for synchrotron emission 
(cf. Rees &; Meszaros 1992, Meszaros & Rees 1993, and references therein). We take the electrons 
to have the energy distribution N{£) oc for fmin < £ < ^'max, i-e., a power law distribution 
with equal energy per decade. This is similar to electron energy spectra inferred from synchrotron 
emission observed in a variety of sources, and the choice of exponent 2 reduces the sensitivity of the 
results to Emm and <?max- Working for the moment in the frame of the shell we adopt the standard 
results for synchrotron emission from a power law distribution of electrons (cf. Pacholczyk 1970, 
§3.4). The spectrum takes the form I,j oc z/^/^(l — e~'^) where r is the optical depth to synchrotron 
self-absorption. For a power law exponent of 2 and frequencies in the range t'c(^^min) ^ ^ ^ 
'^c(^max), we have in the frame comoving with the shell 

r = 0.35 »-(^)„,„,(^)"' cgs (13) 

where w = logiQ(£^inax/^min), Ar is the thickness of the emitting slab, Ue and are the electron 
and the magnetic energy densities, and Vc is the characteristic frequency for synchrotron emission 
of a single electron, given by ydle) ~ 5 x 10^ {£/\mf,c^\) (S/gauss). It follows that the spectrum 
peaks at the frequency i'rn,co for which r = 0.35, rising as v^f^ at lower frequencies and falling off 
as v~^-^ at higher ones. The peak specific intensity is 

hmco = 4.3 X 10-^ ( (^^\ \ ^-^ 1 . (14) 

\gauss/ \5GHz/ Lcm^sHz SteradianJ 

We do not know the magnetic energy density Ug = B'^/Stt or the electron energy density Ue, 
so we express them in terms of dimensionless parameters, 

Ub = Cs^i, Ue = ieUi, (15) 

and conservatively guess that .^e ~ 0.01, well below the equipartition value - the maximum plausible 
range is {rrie/mp) <, <, 0.5; given enough time the magnetic field tends to approach equipartition 
but it is not likely to be that strong from the beginning of the sweeping up phase, so it seems 
reasonable to take .^^ ~ 0.01 as a representative value. We can constrain the range over which the 
electron distribution obeys N{£) oc S"^ by assuming that inverse Compton losses enforce Urad <, . 
This gives an upper bound on w typically in the range 2 to 8, with lower allowing higher w. 

The following is a rather tedious task. First we have to express r, z^co,m and Iu,m,co from eqs. 
(13) and (14) in terms of the fireball parameters given in eqs. (3-5). Next, we have to transform 
to the frame of the Earth-bound observer. The relativistic transformations of the frequency and 
the specific intensity are 

u^ = {i + p)rueo ~ 2ri/eo , I.,® = (1 + Pfr^h,co ~ sr^u^co ■ (i6) 

To go from the specific intensity to a flux density we multiply by the apparent size of the source, 

Fu,®,m ~ ^2^2 ^i^,co,m ■ (17) 

The relation between / and is given in eq. (10). We will restrict our attention to the regime 
ro ' ^ / ^ To, so that te ~ {ft)l{^^o)- 
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After some algebra we find that the time from the beginning of the fireball to the observed 
peak in the flux at frequency is 

1 /2 

and that the corresponding peak flux is 

-2 / F \ 7/8 ^ „ V 1 /S / V 5/8 



F,e,. « 56mJy (iF^)'^'(^)^ 



(19a) 



where 



are dimensionless coefficients of order unity, S is the burst fluence as measured at Earth (E^ = 
AttcIqS), and ^-y = E^/Eq is the fraction of fireball's energy radiated in gamma-rays in in the range 
20 — 2, 000 keV. The bulk Lorentz factor at the time of peak emission at is 

l™*o~3G, ^j^j (^^^j . (21) 

By combining equations (18) and (19) with the spectrum we find that at a fixed frequency the 
fiux rises as prior to t©,m and declines as t^'^^^ thereafter, so that the total time above 
half maximum is about 2t0^m- We can also use (18) and (19) to calculate the peak brightness 
temperature Tf,^rn = c^^i^,co,m/{2i^co,m^)- The result is approximately 10^^ K and is extremely 
insensitive to all parameters of the model. 

3. DISCUSSION 



Observational upper limits to the radio emission from GRBs are available in the literature for 
timescales ^ 10 hours (Baird et al 1975; Cortiglioni et al 1981) and for ^ 5 years (Schaefer 
et al 1989). Baird ct al (1975) used ~ full sky coverage at frequencies Vq < 151MHz to watch 
for strong radio bursts {F^, ^ 100 kJy, duration < 100 s) and later checked for events that had 
occurred within 10 hours of recorded GRBs. Cortiglioni et al (1981) used a similar method at 151 
and 408MHz, attaining a somewhat higher sensitivity (Fj^ ^ 1 kJy) and looking for closer time 
coincidence (< 10 minutes) with GRBs. Schaefer et al (1989) used the VLA at 15, 5, and 1.4GHz 
to conduct a deep search for quiescent emission in IPN GRB error boxes, with sensitivity from 
0.1 to O.SmJy. Other potentially relevant searches have been done by Vaughan and Large (1987) 
and by Amy et al (1989). However, none of these studies found anything above a normal rate 
of background events. Other efforts are underway, but there have been no confirmed detections 
of radio emission from GRBs, nor are we aware of published limits on a timescale of a few days. 
Hanlon, Bennett, and Spoelstra reported a candidate detection (76mJy at 0.6GHz, < ImJy at 
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5GHz) three days after GRB 930309 (lAU Circular 5749; cf. also lAU Circulars 5750, 5755, 5763, 
and 5764 [1993]) but IPN data showed the radio source and GRB positions to be inconsistent. 

None of these experiments provides serious restrictions on our scenario. If gamma-ray bursts 
are at cosmological distances then we expect the peak of the radio emission to follow the burst 
within a week or so, and for the strongest bursts to be up to 0.1 Jy at frequencies of a few GHz. 
It is clear that a dedicated search with the VLA or other sensitive radio instruments should follow 
accurate position determinations from the Interplanetary Network as soon as possible. It would 
be best to observe the rise as well as the fall of the radio transient as this would provide very good 
diagnostics for the burst environment, and in particular for the distance scale. If the bursts arc in 
the galactic halo, their energies are likely to be ~ 10*^^ ergs rather than 10^'^ ergs, and so the time 
scale of the transient would be ~ 10 seconds, i.e. practically simultaneous with the gamma-ray 
bursts (cf. eq. 18). The peak radio flux depends weakly on the distance scale (cf. eq. 19), still the 
galactic halo radio transients might be as strong as ~ 1 Jy. 

It should be noted that there is one major uncertainty in our estimate: the gamma-ray burst 
may be strongly beamed, perhaps to a solid angle as small as ttF^ ~ 10~^. In this case the 
energy required to power the burst may be reduced to a value as small as Eq ~ lO^^erg without 
afi!'ecting the observable properties of the gamma-ray burst, but strongly reducing the radio power 
which peaks at Tradio ~ 3 (cf. eq. 21). The modest amount of the radio beaming may bring the 
equivalent fireball energy up to £"0 ~ lO'*^ erg, but our estimate of the peak radio power is reduced 
by a factor ~ 10* making the radio transient undetectable (cf. eq. 19). To maximize the likelihood 
of detecting a radio transient one should observe at the highest frequency possible, as the total 
radio power increases with frequency (cf. eq. 19) and the radio beaming factor increases as T'^^dio 
(cf. eq. 21). Of course, the higher the radio frequency the more rapid the radio transient (cf. eq. 
18), and the more rapidly must the VLA observation follow the gamma-ray burst. 
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